Systems with superabsorbing states 
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We report on some extensive analysis of a recently proposed model [A. Lipowski Phys. Rev. 
E 60, 6255 (1999)] with infinitely many absorbing states. By performing extensive Monte Carlo 
simulations we have determined critical exponents and show strong evidence that this model is not in 
the directed percolation universality class. The conjecture that this two-dimensional model exhibits 
a dimensional reduction (behaving as one-dimensional directed percolation) is firmly disproven. The 
reason for the model not exhibiting standard directed percolation scaling behavior is traced back 
to the existence of what we call superabsorbing sites, i.e. absorbing sites that cannot be directly 
activated by the presence of neighboring activity in one or more than one directions. Supporting 
this claim we present two strong evidences: (i) in one dimension, where superabsorbing sites do 
not appear at the critical point, the system behaves as directed percolation, and (ii) in a modified 
two-dimensional variation of the model, defined on a honeycomb lattice, for which superabsorbing 
sites are very rarely observed, directed percolation behavior is recovered. Finally, a parallel updating 
version of the model exhibiting a nonequilibrium first-order transition is also reported. 

PACS numbers: 05.10 Ln, 05.50 +q, 05.70 Ph, 05.70 Ln 



I. INTRODUCTION 

Phase transitions separating active from fluctuation- 
free absorbing phases appear in a vast group of physical 
phenomena and models as, for instance, directed perco- 
lation , catalytic reactions Q , the pining of surfaces 
by disorder Q, the contact process Q, damage spread- 
ing transitions 0, nonequilibrium wetting or sand- 
piles HI]. See If and |] for recent reviews. Classifying 
these transitions into universality classes is a first priority 
theoretical task. As conjectured by Janssen and Grass- 
berger some time ago and corroborated by a huge 
number of theoretical studies and computer simulations, 
systems exhibiting a continuous transition into a unique 
absorbing state with no extra symmetry or conservation 
law belong to one and the same universality class, namely 
that of directed percolation (DP). At a field theoretical 
level this class is represented by the Reggeon field theory 
(RFT) 0. ^ ^ 

This universality conjecture has been extended to in- 
clude multicomponent systems a nd systems with in- 
finitely many absorbing states ^3||l4]. On the other 
hand some other, less broad, universality classes of sys- 
tems with absorbing states have been identified in recent 
years. They all include some extra symmetry or conser- 
vation law, foreign to the DP class. For example, if two 
symmetric absorbing states exist (which in many cases 
is equivalent to having activity parity-conservation ) , 
the universality class is other than DP, and the corre- 
sponding field theory differs from RFT [^. A second 
example is constituted by systems with absorbing states 
in which fluctuations occur only at the interfaces sepa- 
rating active from absorbing regions, but not in the bulk 
of compact active regions (examples of this are the voter 



model or compact directed percolation jr^). In this case 
the exponents are also non DP. A third and last example 
is that of systems with many absorbing states in which 
the activity field is coupled to an extra conserved field. 
This type of situation appears, for example, in conserved 
sandpile models, and has been recently shown to define 
a new universality class [p|p^. Apart from these and 
some few other well known examples JT9| , systems with 
absorbing states belong generically into the DP univer- 
sality class. 

Recently, Lipowski has proposed a very simple, bio- 
logically motivated model, exhibiting a continuous tran- 
sition into an absorbing phase, and claimed that this 
model shows a sort of "superuniversality" , i.e. in both 
one and two dimensions the model has the same critical 
exponents, namely those of one-dimensional DP. Con- 
sequently, the system has been hypothesized to show a 
rather strange "dimensional reduction" [ pO[ in two di- 
mensions. This conclusion, if confirmed, would break the 
Janssen-Grassberger conjecture, since it is not clear that 
any new symmetry or extra conservation law is present 
in this model. In what follows we show what are the 
physical reasons why this model does not show directed 
percolation behavior: the presence of what we called su- 
perabsorbing sites is at the basis of this anomalous be- 
havior. We will discuss also how DP can be restored by 
changing the geometry of the lattice on which the model 
is defined. 



II. THE MODEL 

The Lipowski model is defined operationally as follows: 
let's consider a square d-dimensional lattice. At a bond 
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linking neighboring sites, a random variable w = Wij 
is assigned. Different bonds are uncorrelated, and w is 
distributed homogeneously in the interval [0, 1]. At each 
site i one defines as the sum of the four bonds connect- 
ing this site to its four nearest neighbors. If is larger 
that a certain threshold, r (that acts as a control param- 
eter) the site is declared active, otherwise the site is inac- 
tive or absorbing. Active sites are considered unstable; at 
each step one of them is chosen randomly and its four as- 
sociated Wij bond variables are replaced by four freshly 
chosen independent random values (extracted from the 
same homogeneous probability distribution), and time is 
incremented by an amount At = l/n(t), where n{t) is the 
number of active sites at that time. Critical exponents 
are defined as it is customary in the realm of absorbing 
phase transitions 

It is clear that for small values of r, for instance r = 0, 
the system will always be active, while for large enough 
values of r an absorbing configuration (with rt < r for 
all sites i) will be eventually reached. Separating these 
two regimes we observe a critical value of r, Tc, signaling 
the presence of a continuous phase transition. In c? = 1 
rc « 0.4409 H, while for d = 2 we find n = 1.38643(3). 
As bond variables are continuous it is obvious that there 
is a continuous degeneracy of the absorbing state (i.e. 
infinitely many absorbing configurations). 

In the one dimensional case, all the measured critical 
exponents take the expected DP values |23), compatible 
with theoretical predictions for systems with many ab- 
sorbing states |l^,Q. The only discrepancy comes from 
the fact that the spreading exponents 77 and S (see sec- 
tion IIIB for definitions) appear to be non-universal, but 
the combination r] + S coincides with the DP expecta- 
tion. This non-universality in the spreading is however 
generic of one-dimensional swtems with an infinite num- 
ber of absorbing states |25y24| ], and therefore it does not 
invalidate the conclusion that the system behaves as DP. 

In two dimensions the only measured critical exponent 
in [ p6| is the order parameter one, (3, which has been 
reported to take a value surprisingly close to the one di- 
mensional DP expectation, /3 « 0.27 |^^. Based on this 
observation Lipowski claimed that the system exhibits 
a sort of dimensional reduction. This possibility would 
be very interesting from a theoretical point of view and 
elucidating it constitutes the main original motivation of 
what follows. 

Finally let us mention that for spreading experi- 
ments it was found that, as happens generically in two- 
dimensional systems with many absorbing states p^ , p7t , 
the critical point is shifted, and its location depends on 
the nature of the absorbing environment the initial seed 
spreads in. In particular, the annular type of growth de- 
scribed in 1^ in the case of spreading into a favorable 
media is typical of spreading in two-dimensional systems 
with many absorbing states, and it is well known to be 
described by dynamical percolation ||2j,|2^ . 



III. MODEL ANALYSIS 

In order to obtain reliable estimations for /3 and de- 
termine other exponents, we have performed extensive 
Monte Carlo simulations in d = 2 combined with finite 
size scaling analysis, as well as properly defined spreading 
experiments. 

A. Finite size scaling analysis 

We have considered square lattice with linear dimen- 
sion L ranging from 32 to 256. Averages are performed 
over a number of independent runs ranging from 10^ 
to 10^ depending on the distance to the critical point 
and on system size. The first magnitude we measure 
is the averaged density of active sites, p(L,r, t), which 
for asymptotically large times converges to a stationary 
value p{L,r). Observe that for small system sizes the 
system always reaches an absorbing configuration in fi- 
nite time and therefore the only truly stationary state 
is p = 0. In order to extrapolate the right asymptotic 
behavior in the active phase one has to determine p(L, r) 
averaged over the runs which have not reach an absorbing 
configuration. A peculiarity of this system is that its con- 
vergence towards a well defined stationary state is very 
slow, fiuctuations around mean values are extremely per- 
sistent and, therefore, a huge number of runs is needed 
in order to obtain smooth evolution curves. Owing to 
this fact, numerical studies are rather costly from a com- 
putational point of view. The reasons underlying such 
anomalously long lived fluctuations will be discussed in 
forthcoming sections. The maximum times considered 
are 8 x 10^ Monte Carlo steps per spin; this is one order 
of magnitude larger than simulations presented in [ p6[ . 
Near the critical point the relaxation times are very large 
(larger than 10^) and, in order to compute stationary av- 
erages, transient effects have been cut off. We observe the 
presence of a continuous phase transition separating the 
active from the absorbing phase at a value of r « 1.38. 

Assuming that finite size scaling holds (28| in the vicin- 
ity of the critical point point we expect for values of 
r < Tc (i.e. in the active phase) 

p{L,r)r.L-^^''-g{L/A-'-) (1) 

where A = \r — rd- Right at the critical point, this cor- 
responds to a straight line in a double logarithmic plot 
of p{L, r) vs. L. 

In figure 1 it can be seen that, in fact, we observe a 
straight line as a function of log{L) for r = 1.38643(3) 
which constitutes our best estimation of rc- This fi- 
nite size analysis allows us to determine rc with much 
better precision than in the previous estimations [p6| . 
From the slope of the previous log-log plot we measure 
P/i'± — 0.57(2) which is quite far from both, the one- 
dimensional DP exponent P/v± = 0.2520(1), and the 
two-dimensional value 0.795(5). 



2 



50 ^00 300 

FIG. 1. Density of active sites as a function of L (the lin- 
ear system size) for different values of r: from top to bot- 
tom, 1.38630, 1.38640, 1.38643, 1.38645, and 1.38650 respec- 
tively. The straight solid line corresponds to the critical point 
Tc = 1.38643(3). 

We have considered the larger available system size, 
L = 256, and studied the time decay of a fully ac- 
tive initial state for values of r close to rc in the ac- 
tive phase (see figure 2). The stationary values for large 
values of t should scale as p{L,r) ~ A(L)'^. From the 
best fit of our data (see figure 3) we determine both 
rc{L = 256) « 1.38645 and f3 = 0.40(2). 

At the critical point, p{r = rc,t) ~ t^^ . From the 
asymptotic slope of the curve for rc{L — 256) in figure 2, 
we measure 9 — 0.275(15). In this way, we have already 
determined three independent exponents. ^From these, 
using scaling laws, we can determine others, as for exam- 
ple i/± = /3/(/3/j^j_) = 0.69(9) (to be compared with the 
DP prediction 1.09 in ci = 1 and 0.733 in two dimensions 

To further verify the consistency of our results we have 
considered p{L,r) computed for different values of r and 
L, and assumed that p{L,r)L^/^^ , depends on r and L 
through the combination i-'^/'^^A [Q. In figure 4, we 
show the corresponding data collapse which is rather 
good when the previously reported values of fi and v±_ 
are used. In the inset we verify that the data points 
are broadly scattered when one-dimcnsional DP expo- 
nent values are considered, showing that the dimensional 
reduction hypothesis is not valid. Data collapse is nei- 
ther observed using two-dimensional DP exponents; this 
provides a strong evidence that we are in the presence 
of anomalous (non-DP) scaling behavior. Finally, let us 
remark that the observed scaling does not extend over 
many decades for any of the computed steady state mag- 
nitudes. Much better scaling is observed for spreading 
exponents as will be shown in the following section. 
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FIG. 2. Time evolution of the density of active sites for 
L = 256 and different values of r in the active phase, namely, 
from top to bottom r = 1.38143, 1.38402, 1.38527, 1.38587, 
1.38616, 1.38630, 1.38637, and 1.38640 respectively. From the 
slope of the straight dashed line we estimate 6 = 0.275(15). 



B. Spreading experiments, and superabsorbing 
states 

In order to further verify and support our previous con- 
clusion we have performed also spreading experiments as 
it is customarily done in systems with absorbing states 
pO|,P. These consist in locating a seed of activity at 
the center of an otherwise absorbing configuration, and 
studying how it spreads on average in that medium [Q. 
In the absorbing phase the seed dies exponentially fast, 
propagates indefinitely in the active phase, while the crit- 
ical point corresponds to a marginal (power law) propa- 
gation regime 

As said before, it is well established that two- 
dimensional systems with infinitely many absorbing 
states show some peculiarities in studies of the spread- 
ing of a localized activity seed. The absorbing environ- 
ment surrounding the seed can either favor or un-favor 
the propagation of activity depending on its nature (see 
p^ , p7[ and references therein). For the, so called, nat- 
ural initial conditions |0| the critical point for spreading 
coincides with the bulk critical point, and standard DP 
exponents are expected. In order to generate such natu- 
ral configurations one could start the system with some 
highly active configuration and run the system right at 
the critical point; once it reaches an absorbing configu- 
ration it can be taken as a natural or self-generated envi- 
ronment for spreading. An alternative, more efficient way 
of proceeding, inspired in sandpile systems is as fol- 
lows. One considers an arbitrary absorbing configuration 
and runs a spreading experiment. Once the epidemic (or 
"avalanche" in the language of self-organized criticality 
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A(L=256)=(r^XL=256)-r) 
FIG. 3. Stationary density of active sites as a function of 
the distance to the critical point, for L — 256 and different val- 
ues of r in the active phase (the same values reported in figure 
2). The best fit gives 13 = 0.40(2) and rc{L = 256) « 1.38645. 
Filled (empty) circles are used to represent scaling (not scal- 
ing) points. 

is over, one considers the newly reached absorb- 
ing configuration as initial state for a new spreading ex- 
periment avalanche. After iterating this process a num- 
ber of times the system reaches a statistically stationary 
absorbing state: the natural one (see |^ and references 
therein) . Using this absorbing state for spreading leads to 
DP exponent values (and critical point) in systems with 
many absorbing states as for example the pair contact 
process |25 3^. 

By following this procedure we have found a very pe- 
culiar property of this model, that we believe to be at 
the basis of its deviating from DP. If the initial seed is 
located for all avalanches in the same site (or small group 
of localized sites), as is usually the case, after a relatively 
small number of avalanches the system reaches an absorb- 
ing configuration such that it is impossible to propagate 
activity for any possible forthcoming avalanche beyond 
a certain closed contour. For example, configurations 
as the one showed in figure 5a are generated. The four 
sites at the center are the ones at which activity seeds 
are placed in order to start avalanches. White sites are 
active and grey ones are absorbing. At each marked-in- 
black site, the sum of the three (black) bonds connect- 
ing it to sites other than a central one is smaller than 
Tc — 1 = 0.38643(3). In this way, regardless the value 
of the bond connecting the site to the central region the 
site remains inactive: it is a superahs orbing site. The 
existence of "inactive forever" sites have been already 
pointed out by Lipowski ||2^,Q. In the configuration 
showed in figure 5a activity cannot propagate out of the 
"fence" of superabsorbing sites: the cluster of superab- 
sorbing sites will remain frozen indefinitely, and activity 



FIG. 4. Data collapse for the density of active sites: 
p'(L,A) = p(L, A)L'^/''^ and A' = AL^/^-l. Using the ob- 
tained exponent values, l3/v± « 0.57, and u±_ ~ 0.69, a rea- 
sonably good data collapse is observed. In the inset we show 
an attempt to collapse data using one- dimensional DP expo- 
nent values. There is no evidence of scaling neither in this 
case nor using two-dimensional DP exponents. 

cannot possibly spread out. All avalanches will neces- 
sarily die after a few time steps. This type of blocking 
structure is quite generic, and appears in all experiments 
after some relatively short transient. 

In conclusion, this way of iterating spreading experi- 
ments leads always to blocking closed configurations of 
superabsorbing sites instead of driving the system to a 
natural absorbing configuration. 

Observe that some activity put out of a blocking fence 
of sites in figure 5a could well affect any of the exter- 
nal bonds of the superabsorbing sites (the dangling black 
bonds in figure 5a), converting the corresponding site to 
an absorbing or even an active one. Therefore, in order 
to overcome this difficulty of the frozen blocking configu- 
rations and be able to perform spreading experiments in 
some meaningful way, we iterate avalanches by locating 
the initial seed at randomly chosen sites in the lattice. 
In this way there is always a nonvanishing probability 
of destroying blocking "fences" by breaking them from 
outside as previously discussed. Measurements of the 
different relevant magnitudes are stopped when the sys- 
tem falls into an absorbing configuration or alternatively 
whenever a linear distance L/2 from the avalanche origin 
is reached. Observe that in the second case the dynamics 
has to be run farther in order to reach a new absorbing 
configuration at which launching the next avalanche. 

We monitor the following magnitudes: the total num- 
ber of active sites in all the runs as a function of time N{t) 
(we also estimate Ng [t) defined as the average number of 
active sites restricted to surviving runs), the surviving 
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FIG. 5. Different frozen configurations of superabsorbing 
(black) sites. White (grey) color stands for active (absorb- 
ing) sites, (a) Blocking configuration for spreading from the 
central cluster of four sites. Black sites cannot change their 
state whatever the state or dynamics inside the cluster might 
be. Black bonds remain also frozen, (b) Spanning frozen 
cluster of superabsorbing sites, (c) Almost-frozen cluster of 
superabsorbing sites. This, and analogous structures, can be 
destabilized from the outside corners. 

probability P{t), and the average square distance from 
the origin, (t) . At the critical point these are expected 
to scale as N{t) - t'', P{t) ^ and R'^it) Resuhs 
for this type of measurements are reported in figure 6. 

We obtain rather good algebraic behaviors at the pre- 
viously estimated critical point, rc, confirming that the 
iteration-of-avalanches procedure leads the system to a 
natural absorbing environment. Slightly subcritical (su- 
percritical) values of r generate downward (upward) cur- 
vatures in this plot for all the four magnitudes. Our best 
estimates for the exponents at criticality are: z = 0.96(1), 
rj = 0.05(1), S = 0.66(1) (see table 1). To double check 
our results we also plot Ns{t), which is expected to scale 
with an exponent i] + S. An independent measurement 
of its slope in the log-log plot gives t] + S — 0.71(1), in 
perfect agreement with the previously obtained results. 

We can use these values to verify the hyperscaling re- 
lation MM 
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viously determined value = 0.275(15). 

One more check of the consistency of our result s by 
using scaling laws is the following. As z = 2i^^/i/|| p9[ , 
we can estimate i/y from z and i^^. Then, using and 
the fact that — l3/iy\\ we obtain 6 ~ 0.27(1), again in 
excellent agreement with the directly measured value. 

In table 1, we present the collection of exponents and 
compare them with DP values in both one and two di- 
mensions 1 29 . There is no trace of dimensional reduction: 



this model does not behave, at least up to the scales we 
have analyzed, as any other known universality class. 



C. More about superabsorbing states 

Let us recall our definition of superabsorbing states. 
A site, three of whose associated bonds take values such 
that the sum of them is smaller that r — I, cannot be 
activated from the remaining direction by neighboring 
activity. We say that this site is superabsorbing in that 
direction (or it is in a superabsorbing state). A site can 
be superabsorbing in one or more than one directions. 
Still a site in a superabsorbing state can obviously be ac- 
tivated by neighboring activity in any of the remaining 
directions (if any). 

Having stated the existence of frozen clusters in stan- 
dard spreading experiments (when initialized from a fixed 
localized set of sites), one may wonder whether there 
are similar frozen structures in simulations started with 
an homogeneous initial distribution of activity, or in the 
modified type of spreading experiments we have just used 
(i.e. allowing the initial seed to land at a randomly cho- 
sen site) in the neighborhood of the critical point. 

In principle, for any finite lattice, the answer to that 
question is affirmative. In figure 5b we show the shape of 
a frozen cluster of superabsorbing sites: any of the sites 
in it is superabsorbing with respect to the corresponding 
outward direction, and it cannot be "infected" from any 
of the other directions as neighboring sites are similarly 
superabsorbing. If a cluster like that is formed (or put by 
hand on the initial state) it will remain superabsorbing 
forever. However the probability to form such a perfectly 
regular chain is extremely small for large system sizes. 
Observe also that in order to have a completely frozen 
two-site broad band structure it has to be unlimitedly 
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0.275(15) 


0.05(1) 
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DP, d = 1 


0.276 
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0.583 


0.795 


0.450 


0.229 


0.450 


1.132 



Substituting the found values for z and -t- (5 we obtain 
d K. 0.25(2), compatible within error bars with the pre- 



TABLE I. Exponent values for the two dimensional 
Lipowski model and directed percolation in both one and two 
dimensions. Figures in parenthesis denote statistical uncer- 
tainty (note that error-bars are statistical errors coming from 
power-law fittings, and therefore do not include eventual sys- 
tematic corrections to scaling). 
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FIG. 6. Numerical results for spreading experiments. R^{t) 
(topmost curve), Ns{t) (second curve from above), N{t) 
(third curve from above), and P{t) (bottom curve). ^From 
the slopes we estimate z — 0.96(1) and rj + 5 = 0.71(1), 
r; — 0.05(1) and 5 = 0.66(1) respectively. 




FIG. 7. Snapshot of a configuration in a 32 x 32 lattice in 
the stationary regime for a value of r close to the critical point. 
White color denotes activity, black corresponds to superab- 
sorbing sites, while grey stands for absorbing sites. Observe 
that superabsorbing sites percolate through the lattice. 



long (or closed if periodic boundary conditions are em- 
ployed). If instead it was finite, then sites at the corners 
would be linked to two external susceptible-to-change 
bonds and, therefore, themselves would be susceptible 
to become active: they would not be blocked forever. In 
this way any finite structure of superabsorbing sites in 
the square lattice is unstable: it can be eaten up (though 
very slowly) by the dynamics, and is therefore not fully 
frozen. For instance, the cluster of superabsorbing sites 
represented in figure 5c is almost-frozen, but not really 
frozen as it may lose its superabsorbing character from 
the outside corners as previously described. Analogously, 
any other cluster shape of superabsorbing sites may be 
destabihzed from its outside corners. 

In conclusion, frozen clusters of superabsorbing sites 
do not appear spontaneously. However, almost-frozen 
regions do appear and may have extremely long life 
spans, specially close to the critical point where activ- 
ity is scarce, and therefore the possibility of destabilizing 
them is small. In order to give an idea of how frequently 
superabsorbing sites appear we present in figure 7 a snap- 
shot of a typical system-state near the critical point. 

White corresponds to active sites, while the remaining 
sites are absorbing: in black we represent superabsorbing 
(in one or more than one directions) sites, while simple 
absorbing (not-superabsorbing) sites are marked in grey 
color. Observe that superabsorbing sites are ubiquitous; 
in fact they percolate through the system. Among them, 
about one forth are superabsorbing in all four directions. 

Even though none of the clusters of superabsorbing 
sites is completely frozen, and in principle, activity could 
reach any lattice site, the dynamics is glassy ]34| in some 
sense. For instance, imagine an active region separated 



from an absorbing region by a line of superabsorbing- 
in-thc-dircction-of-the- activity sites. In order to reach 
the absorbing region, activity has to circumvent the su- 
perabsorbing barrier. But near the critical point, where 
activity is scarce, barriers of superabsorbing sites are in- 
tertwisted among them forming structures that, even if 
not completely frozen, are very unlikely to be infected: 
activity has to overcome them progressively in order to 
reach the interior of superabsorbing regions. This re- 
sembles some aspects of glassy systems for which degrees 
of freedom are hierarchically coupled and, at observable 
timescales, they may appear effectively frozen p^ . 

This phenomenology is certainly very different from 
DP, and it is the reason why the relaxation towards sta- 
tionary states is so slow, and why deviations from mean 
values are so persistent in numerical simulations. In par- 
ticular, as superabsorbing regions are long lived, the time 
required for the system to self-average is very large, and 
as near the critical point the probability of reaching an 
absorbing state is large, in practice, the system does 
not have the time to self-average. Consequently, a huge 
amount of independent initial states and runs have to 
be considered in order to measure smooth well behaved 
physical magnitudes |3^. We strongly believe that this 
type of pathological dynamics is responsible for the de- 
parture of the Lipowski model from the DP universality 
class in two dimensions. 

At this point one might wonder whether the one- 
dimensional version of this model is essentially different. 
Or in other words, why (one-dimensional) DP exponents 
are observed in d = 1 |^^?. The answer to this question 
is not difficult if one argues in terms of superabsorbing 
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sites. First of all notice that in d = 2, Tc > 1. This means 
that just by changing one bond, whatever the value of the 
output is, the site can stay below threshold if the other 
three bonds sum less than rc — 1; this is to say superab- 
sorbing states do exist at criticality. However in d — 1, 
Tc = 0.4409 < 1. In this case by changing one bond 
value it is always possible to activate the corresponding 
site: superabsorbing sites do not exist in d = 1 at the 
critical point |^^. Once the "disturbing" ingredient is 
removed from the model, we are back to the DP class as 
general principles dictate. 



D. The honeycomb lattice 

In order to further test our statement that superab- 
sorbing states are responsible for the anomalous scaling 
of the two-dimensional Lipowski model, we have stud- 
ied the following variation of it. We have considered the 
model defined on a honeycomb lattice (with three bonds 
per site), and performed Monte Carlo simulations. In this 
case there is the (geometrical) possibility of having com- 
pletely frozen clusters of superabsorbing sites (see figure 

. . . 

The main geometrical difference from the previous case 

comes from the fact that here cluster-corners are linked 
only to one external bond, and therefore are more prone 
to form frozen clusters. In principle, before performing 
any numerical analysis, there are two alternative pos- 
sibilities: either the critical point is located at a value 
of r smaller than 1 or larger than 1. In the first case, 
there would be no superabsorbing site (in analogy with 
the one-dimensional case); in the second case pathologies 
associated with superabsorbing sites should be observed. 
The case rc — 1 would be marginal. Finite size scaling 
analysis indicate the presence of a continuous phase tran- 
sition located at r w 1.0092 (very nearby the marginal 
case, but significatively larger than r — I). 

For Monte Carlo simulations, we have employed lat- 
tices of up to a maximum of 256 x 256 sites. All the ob- 
served phenomenology is perfectly compatible with two- 
dimensional DP behavior. The dynamics does not show 
any of the anomalies described for the square lattice case. 
In particular, from the dependence of the stationary ac- 
tivity density on system size we evaluate (3/i' — 0.80(1); 
from the time decay at criticality 9 — 0.45(1), and fi- 
nally P — 0.57(2); fully confirming consistency with two- 
dimensional DP behavior. This result seems to be in 
contradiction with the two alternative possibilities pre- 
sented above. Let us now discuss why this is the case. 

As the coordination number is 3 in this case, the sum 
of two bond- values has to be smaller than Tc — 1 w 0.0092 
in order to have a superabsorbing site in the direction of 
the remaining bond at criticality. As the two bonds are 
independent random variables, the probability of creat- 
ing a superabsorbing site if the two of them are changed, 
is fewer than 0.5%, and the probability to generate frozen 



clusters (composed by six neighboring superabsorbing 
sites as shown in figure 8), is negligible at the critical 
point. In fact, we have not been able to observe any of 
them in our simulations. This means that one should 
study extremely large system sizes and extraordinarily 
long simulations in order to see anomalies associated with 
superabsorbing sites, otherwise, for any feasible simula- 
tion the behavior is expected to be DP-like. The obser- 
vation of DP exponents in this case strongly supports 
the hypothesis that superabsorbing states are at the ba- 
sis of the anomalous behavior of the model on the square 
lattice. 

However, strictly speaking, the system should exhibit 
a (unobservable) first-order phase transition at r = 1 in 
the thermodynamic limit. Indeed, for values of r larger 
than 1 there is a finite, though extremely small, proba- 
bility of creating frozen clusters of superabsorbing sites 
(as the one in figure 8). As this is an irreversible process, 
after some (divergently long) transient there would be a 
percolating network of frozen clusters of superabsorbing 
sites, and the only possible stationary state would be an 
absorbing one with zero activity. On the other hand, for 
values of r smaller than unity, the probability of creat- 
ing superabsorbing sites is strictly zero, and there will 
be a nonvanishing density of activity. As the density at 
r = 1, almost independent of system size, is p « 0.18, the 
transition is expected to be discontinuous, and therefore 
the DP transition observed in our simulations is merely 
a finite size effect, and should disappear for large enough 
sizes and long times. In any case, this first order transi- 
tion is unobservable computationally. 



IV. SUMMARY 

Summing up, we have shown that the two-dimensional 
Lipowski model does not belong to any known universal- 
ity class. We have measured different critical exponents 
by running Monte Carlo simulations started from homo- 
geneous initial states and also by performing spreading 
experiments. In any case, we find absolutely no trace of 
dimensional reduction, and there is neither evidence for 
the system to behave as two-dimensional DP. Instead, a 
novel scaling behavior is observed. The main relevant 
physical ingredient of this class is the presence of super- 
absorbing sites, and almost-frozen clusters of superab- 
sorbing sites which slow down enormously the dynamics. 

The previous conclusion is strongly supported by two 
other observations: (i) the regular DP behavior observed 
in the one-dimensional version of the model for which su- 
perabsorbing states do not appear at criticality, and (ii) 
the two-dimensional DP behavior observed for the two- 
dimensional model defined on a honeycomb lattice, for 
which the probability of generating superabsorbing sites 
at criticality is almost negligible. 

In general, superabsorbing sites can either arrange into 
completely frozen clusters or not depending on dimen- 
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FIG. 8. Frozen cluster in the honeycomb lattice. This type 
of frozen structure of superabsorbing sites remains indefinitely 
superabsorbing at the critical point. Black: superabsorbing 
sites. Grey: absorbing sites. White: active sites. 

sionality, coordination number and other system details. 
Let us distinguish three main cases: 

(1) When completely frozen clusters of superabsorbing 
sites appear below (or above) a certain value of the con- 
trol parameter but not above (below), first order transi- 
tions are expected (as occurs in the multiplicative model 
discussed in appendix 2 ||3^). 

(2) If completely frozen clusters do not appear at crit- 
icality, but instead almost-frozen clusters are present, 
we expect anomalous behavior (as occurs in the original 
Lipowski model [2^]). 

(3) If neither frozen nor almost-frozen clusters are ob- 
served at criticality (as is the case for the one dimensional 
version of the model [p3|) we expect standard directed 
percolation behavior. 

Two possible follow-ups of this work are: 

(1) It would be worth studying in more realistic situ- 
ations as, for instance, in surface catalysis (dimer-dimer 
or dimer-trimer) models |13| whether effects similar to 
those described in this paper play any relevant role. In 
particular, for those models depending upon lattice and 
particle geometry there are cases in which activity can- 
not propagate to neighboring regions, but is constrained 
to evolve following certain directions or paths. It would 
be rather interesting to sort out whether anomalies re- 
ported for those models ||l^ are related to the existence 
of superabsorbing states. 

(2) From a more theoretical point of view, an inter- 
esting question is: what is the field theory or Langevin 
equation capturing the previously described phase transi- 
tion with superabsorbing states? and, how does it change 
with respect to Reggeon field theory?. Establishing what 
this theory looks like, would clarify greatly at a field the- 
oretical level the effect of superabsorbing states on phase 



FIG. 9. Order parameter as a function of r in the case of 
parallel updating. The transition appears to be discontinuous 
in this case, exhibiting also a hysteresis loop. 



transitions, and would permit to shed some light on the 
degree of universality of this anomalous phenomenology. 
Our guess is that a Reggeon field theory 00 with 
a spatio-temporal dependent anisotropic Laplacian term 
(which, for example, would enhance, un-favor or forbid 
diffusion from certain sites in certain directions) could 
be a good candidate to describe this new phenomenol- 
ogy. Analogously to what happens in field theoretical 
descriptions of other systems with many a absorbing 
states p^ , p^ , the inhomogeneous Laplacian-term coeffi- 
cient should be described by a second physical field cou- 
pled to the activity field in such a way that its fluctua- 
tions would vanish upon local absence of activity. Fur- 
ther pursuing this line of reasoning is beyond the scope 
of the present paper. As long as this program has not 
been completed, is not safe to conclude unambiguously 
that the anomalies described in this paper are relevant 
in the limit of extremely large times and system sizes. 

APPENDIX I 

As an alternative attempt to speed up the dynam- 
ics, and examine further some properties of the two- 
dimensional model, we have implemented the microscopic 
dynamics replacing the original sequential updating by a 
synchronous or parallel one, i.e. all active sites are "de- 
activated" simultaneously at each Monte Carlo step, and 
all their associated bonds are replaced by new random 
variables simultaneously. In this way, as random num- 
bers do not have to be extracted to sequentially select 
sites, the dynamics is largely accelerated. For this mod- 
ified dynamics, we have examined some relatively large 
system sizes, L — 256, and concluded that the nature of 
the transition is changed with respect to the sequential 
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updating case: in this case the transition is first order and 
critical exponents cannot be defined. To show that this 
is the case, in figure 9 we present the stationary activity 
curve. The upper curve corresponds to simulations per- 
formed taking an initial activity-density equal to unity. 
On the other hand, the lower curve is obtained by start- 
ing the system with a natural absorbing configuration, 
and activating on the top of it a small percentage of sites 
(about a ten percent). 

For values of r in the interval 1.545, w 1.555] the 
system reaches different states depending upon the ini- 
tial condition. The presence of a hysteresis loop is a trait 
of the transition first-order nature. First order absorb- 
ing state transitions have been observed in other contexts 
pTf . However, we caution the reader that, as the transi- 
tion is found to occur at a value of r for which the proba- 
bility of creating superabsorbing sites is very large (much 
larger than in the sequential case), and the dynamics is 
therefore extremely anomalous and slow, it could be the 
case that the first order character of the transition is only 
apparent. Extracting clean, conclusive results in the crit- 
ical zone is a computationally very expensive task, that 
we have not pursued. 

APPENDIX II 

Very recently, Lipowski has introduced a multiplica- 
tive version of his model on the square lattice in which 
sites are declared active if the product of the four adja- 
cent bonds is smaller than a certain value of the control 
parameter r Bonds take uncorrelated values in the 
interval [—0.5,0.5] extracted from a homogeneous distri- 
bution. For values of r smaller than r = there is a finite 
(not small) probability to generate superabsorbing sites. 
In this case, it is not difficult to see that isolated super- 
absorbing sites remain frozen forever. In analogy with 
the discussion of the honeycomb-lattice model, a first or- 
der transition is expected at rc = (as discussed also in 
P2[). However, in this the probability to cre- 

ate superabsorbing sites is not negligible, the first order 
transition is actually observable. Based on a numeri- 
cal measurement of /3 Lipowski concludes that the model 
shares first-order properties with second-order features. 
In particular, the transition is clearly shown to be dis- 
continuous, there is no diverging correlation length, but 
/3 is claimed to be however in the two-dimensional DP 
class. Our guess is that this apparent puzzle is simply 
due to a numerical coincidence, and that in fact there 
is no trait of any second-order phase transition feature 
(observe that the fit for beta in spans for less than 
half a decade in the abscise of the log- log plot). 
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